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ABSTRACT 

Context. Low metallicities, large gas-to-star mass ratios, and blue colors of most low surface brightness (LSB) galaxies 

imply that these systems may be younger than their high surface brightness counterparts. 

Aims. We seek to find observational signatures that can help to constrain the age of blue LSB galaxies. 

Methods. We use numerical hydrodynamic modelling to study the long-term (~ 13 Gyr) dynamical and chemical 

evolution of blue LSB galaxies adopting a sporadic scenario for star formation. Our model galaxy consists of a thin gas 

disk, stellar disk with a developed spiral structure and spherical dark matter halo. Our models utilize various rates of 

star formation and different shapes of the initial mass function (IMF). We complement hydrodynamic modelling with 

population synthesis modelling to produce the integrated B — V colors and Ha equivalent widths. 

Results. We find that the mean oxygen abundances, B — V colors, Ha equivalent widths, and the radial fluctuations 

in the oxygen abundance, when considered altogether, can be used to constrain the age of blue LSB galaxies if some 

independent knowledge of the IMF is available. Our modelling strongly suggests the existence of a minimum age for 

blue LSB galaxies. Model B — V colors and mean oxygen abundances set a tentative minimum age at 1.5-3.0 Gyr, 

whereas model Ha equivalent widths suggest a larger value of order 5-6 Gyr. The latter value may decrease somewhat, 

if blue LSB galaxies host IMFs with a truncated upper mass limit. We found no firm evidence that the age of blue LSB 

galaxies is significantly smaller than 13 Gyr. 

Conclusions. The age of blue LSB galaxies may vary between 1.5-6 Gyr and 13 Gyr, depending on the physical conditions 
in the disk that control the form of the IMF and the rate of star formation. A failure to observationally detect large 
radial fluctuations in the oxygen abundance of order 0.5-1.0 dex, which, according to our modelling, are characteristic 
of (l-2)-Gyr-old galaxies, will argue in favour of the more evolved nature of blue LSB galaxies. 

Key words. ISM:galaxies - galaxies: abundances - galaxies: evolution 



1. Introduction 

Low surface brightness (LSB) galaxies are the galaxies 
whose central surface brightness is much fainter than the 
Freeman value /xg = 21.65 ± 0.30 mag arcsec -2 . These 
galaxies are thought to represent a signi ficant fraction of 
the g a laxy number densi t y in the universe (iMcGaugh et all 
119951 IP 'Neil fc Bothunl . l2000t iTrachternach et all 1200(1 . 
LSB galaxies show a wide spread of colours ranging 
from red (B — V ~ 1.7) to very blue (B — V ~ 
0.2) (jO'Neil et all Il997t ). but the most common type 
seems to be blue LSB galaxies: late-type, disk dom- 
inated spirals with central surface brightness /j,b ii 
23 mag arcsec -2 and colors lying in the range (B — 
V) = 0.3 — 0.7. These blue LSB galaxies are among the 
most gas-rich (up to 50% of the total baryonic mass) 
and metal-de ficient (about 5%-20% of the solar metallic- 
ty) galaxies (McGaugh, 1994; Ronnback & Bcrgvall, 



McGaugh fc Blokl. 119971; Ide Blok fc van der HulstT 



Kuzio de Narrav et all 120041 ) 



1995; 



1998; 



Low metallicities and large gas-to-star mass ratios may 
be interpreted in various ways. For instance, blue LSB 



Send offprint requests to: E. I. Vorobyov 

* on leave from Department of Physics, South 
University, and Special Astrophysical Observatory 



galaxies may have formed recently, at redshifts close to 
z = 0. However, this can confront the hierarchical scenario 
of large-scale structure formation, within which low-mass 
galaxies form first at redshifts close to the onset of the reion- 
ization epoch, z ~ 10. Therefore, it is more likely that blue 
LSB galaxies have formed at large z but have had a delayed 
onset of major star formation epoch. In this conflict 
with the hierarchical structure formation is relaxed. Indeed, 
blue LSB galaxies show a wide range of masses with an ex- 
cess of dwarf systems wit h luminosities of order 10 7 — 10 8 L© 
(jSpravberrv et all I1997D . Such dwarf systems are known 
to have a younger mean stellar age than their more mas- 
sive counterparts - a phe nomenon known as 'downsizing' 
(e.g. lThomas et all 120051) . Another possible interpretation 
of the apparent unevolved nature of blue LSBs is that they 
do not form late, nor have a delayed onset of star formation, 
but simply evolve slower than their high surf ace brightness 
counterparts (e.g. Ivan den Hoek et all 120001 ). 

There have been many efforts to infer the age of 
blue LSB galaxies from the properties of stellar popula- 
tion. Broadband photometric studies, complemented by Ha 
emission line data and synthetic stellar population code 
modeling, predict quite a wi de range for the ages o f blue 



Federal LSB ga laxies: from 1-2 Gyr (IZackrisson et all 20051 ) to 7 



9 Gyr (|Padoan et all 119971 : iJimenez* et all 1 1998I L Rather 
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young ages of blue LSB galaxies (< 5 Gyr) follow also 
from their (V — I) colors and HI content (|Schombert et all . 
2001). The comparison of measured spectral energy distri- 
butions in a sample of LSB galaxies with the synthetic spec- 
tra (derived assuming an exponentially declining SFR) by 
iHaberzettl et al.l (|2005T ) suggests the ages of the dominant 
stellar population between 2 and 5 Gyr. 

In this paper, we perform numerical hydrodynamic 
modelling of the dynamical, chemical, and photometric 
evolution of blue LSB galaxies with the purpose to con- 
strain their typical range of ages. The cornerstone sup- 
position of our modelling, motivated by observational evi- 
dence, is a spora dic nature o f star f ormation in blue LSB 
galaxies. Indeed. iMcGaugbl (|l994f) . Ide Blok et all (|l995f) . 
iGerretsen fc de Blok ( 1999f) . and others argue that the pri- 
mary cause for the blue colors is the young age of the dom- 
inant stellar population. This can take the form of young 
local bursts of star formation superimposed (or not) on 
a continuous low-rate star fo rmation. Available Ha im- 
ages of LSB galaxies ( see e.g. McGaugh fc Bothunl . 11994 



time scales for these processes in LSB galaxies are not well 
known and numerical simulations are needed to determine 
their importance. Therefore, the distribution of heavy ele- 
ments in the disk of an LSB galaxy may be characterized 
by radial variations, the typical amplitude of which is ex- 
pected to diminish with time. It is our purpose to determine 
whether these radial abundance fluctuations can be used to 
constrain the ages of LSB galaxies. 

We note that mixing in the azimuthal direction is 
strongly enhanced by differential rotation. Namely, the size 
of a mixing eddy grows with time due to differential ro- 
tation as Ax ~ (u^)' j 1 / 2 ^ 3 / 2 much faster than due t o 
hydrodynamic mixing Ax cx t (jScalo fc Elmegreenl . 1200 J ), 
where (u<f 1 ) r is the radial derivative of the rotation velocity, 
D is the diffusion coefficient. Differential rotation therefore 
homogenizes metals in the azimuthal direction faster than 
in the radial one and the radial variations remain for a 
longer time. 



iMcGaugh etail fl995t lAuld et all . 120061) appearto sup- 3 - Model LSB galaxy 



port this so-called sporadic star formation scenario. Current 
star formation is localized to a handful of compact regions. 
There is little or no diffuse Ha emission coming from the 
rest of the galactic disk. Randomly distributed star forma- 
tion sites, along with inefficient spatial stirring of heavy el- 
ements, give rise to radial fluctuations in the oxygen abun- 
dance on spatial scales of order 1 — 2 kpc. The resulted 
fluctuation spectrum is age-sensitive and can be potentially 
used to set order- of -magnitude constraints on the age of 
blue LSB galaxies. However, the age prediction may be con- 
siderably improved by using the mean oxygen abundances 
and Ha equivalent widths provided that some independent 
knowledge of the initial mass function is available. 

The paper is organized as follows. In the next Section 
we discuss qualitatively the processes that mix heavy ele- 
ments in disks of LSB galaxies. Section [3] details our the- 
oretical model, while Section U describes our observations. 
In Sections [5] and [6] we discuss our numerical results, and 
Section [7] presents main conclusions. 

2. Mixing of heavy elements in LSB disks 

In high surface brightness (HSB) galaxies mixing of heavy 
elements is provided by multiple shock waves from super- 
nova explosions. It requires approximately 100 Myr to ho- 
mogenize heavy elements in the warm neutral interste l- 
lar medium of HSB galaxies (|de Avillez fc MacLowl . l2002D . 
LSB galaxies, on the other hand, are characterized by the 
star formation rate (~ Q.IMq yr _1 ) that is almost two or- 
ders of magnitude lower than the Milky Way value, w 5Mq 
yr _1 . Consequently, the frequency of shock waves in the 
interstellar medium of LSB galaxies and the correspond- 
ing volume filling factor occupied by the shock processed 
gas can be as small as * s/l < (200 Myr) -1 and f sh < 0.004, 
respectively. For comparison, the corr esponding values in 
the M ilky Way are ^ sh ~ (5 Myr)" 1 dDraine fc Salpeterl . 
Il979h and f sh ~ 0.2 (jde Avillez fc BreitschwerdtL 12005( 1 . In 
these conditions, homogenization of heavy elements in LSB 
galaxies may require a considerably longer time scale than 
100 Myr, a typical value for HSB galaxies. Differential ro- 
tation of the galactic disk and radial convection driven by 
the spiral stellar density waves may in principle be efficient 
in mixing of heavy elements in LSB galaxies. However, the 



3.1. Initial configuration 

LSB galaxies show a wide range of morphological types, 
ranging from giant spirals to dwarf irregular galaxies, from 
early-type red galaxies to late-type blue ones. In this pa- 
per we focus on late-ty pe blue LSB galaxies . Acco rding to 
the observational data (Ivan den Hoek et~aL1 . l200(^ . a typi- 
cal blue LSB galaxy has a dynamical mass of order 10 10 Mq 
and HI mass of order 2 x 10 9 Mq. We closely follow these 
estimates when constructing our model galaxy. Our model 
LSB galaxy consists of a thin gas disk, which evolves in the 
external gravitational potential of both the spherical dark 
matter halo and extended stellar disk. LSB galaxies often 
feature a few ill-defined spiral arms. Therefore, we assume 
that our model stellar disk has developed a two-armed spi- 
ral structure. Below, we provide a more detailed explana- 
tion for each of the constituents of our model galaxy. 

The velocity field of most LSB galaxies is better 
fitted by dark matter models described by a modi- 
fied isothermal sphere rath er than by a NFW profile 
(jKuzio de Narav et all 12006( 1 . Therefore, the spherical dark 
matter halo in our model is modelled by a radial density 
profile described by the modified isothermal sphere 



Ph 



PhO 



(1 +r/r h ) 2 



(1) 



where pho = 6.0 x 10~ 3 M pc~ 3 and r h = 5.7 kpc are the 
central volume density and characteristic scale length of the 
dark matter halo, respectively, and r is the radial distance. 
The adop ted values of Pho and are similar to those ob- 
tained bv lKuzio de Narav et all (120061) for UGC 5750. The 
dark matter halo mass in the inner 15 kpc is approximately 
2.0 x 10 10 Mq. The radial gravity force (per unit mass) due 
to the spherical dark matter halo in the plane of the disk 
can be written as 



dr 



= 47rGpho^h [r/fb — arctan 



(r/r h )](^y 



(2) 



The stellar component of our model galaxy consists of 
two parts - an axisymmctric stellar disk and a two-armed 
spiral pattern. The former is described by a radially declin- 
ing exponential profile 



Est = E s0 exp(-r/r s ), 



(3) 



Vorobyov et al.: Age of blue LSB galaxies 



3 



where E s o = 30 M@ pc -2 is the central stellar surface den- 
sity and r a = 4 kpc is the radial scale length of the stellar 
disk . The adopted value f or r s is typical for blue LSB galax- 
ies (|de Blok et all |1995|) . We note that our model stellar 
disk serves as a proxy to the real stellar disk that is expected 
to form during the evolution of our model galaxy. The mass 
of the stellar disk in the inner 15 kpc is M s = 2.6 x 10 9 M Q . 
For numerical purposes, the gravitational potential of the 
axisymmctric thin stellar disk in the plane of the disk is 
calculated as 

$ st (r, z = 0) = -7rGS s0 r [I {v)Ki{v) - h(y)K (y)} , (4) 

where /„ and K n are modified Bessel functions of the first 
and second kinds, respectively, and y = r/(2r 8 ). 

The non-axisymmetric gravitational potential of stellar 
spiral arms is described in polar coordinates (r, </>) by a run - 
ning density wave (see e.g. IVorobvov fc Shchekinovi . [2005) 

$ sp (r, </>) = -C(r) cos [m(cot(i) ln(r/r sp ) + - Q sp t)] , (5) 

where C(r) is the radially varying amplitude, i is the pitch 
angle, r sp is the characteristic radius of the spiral at <j> = 0, 
m = 2 is the number of spiral arms, and Q sp is the angular 
velocity of the spiral pattern. In the following we adopt 
r sp = 6 kpc and il sp = 5.5 km s _1 kpc. The pitch angle is 
set to 25°, typical for the late-type spiral galaxies. 

The amplitude C(r) of the spiral stellar gravitational 
potential <& sp determines the response of the gas and con- 
sequently the appearance of a spiral pattern in the gas 
disk. We adopt the following expression for the amplitude 
C(r) = Co(V) Q ( r *' . Here, Co(r) is a linear function of radius 
which has a value of at r = kpc and equals 1.8 x 10~ 3 
(in dimcnsionless units) at r = 17 kpc. The amplitude is 
set to zero at the galactic center to ensure that the gravi- 
tational potential of the stellar spiral density wave dimin- 
ishes at small galactocentric distances. Otherwise, even a 
small-amplitude spiral gravitational potential would pro- 
duce strong azimuthal gravitational forces near the galactic 
center due to converging radial gridlines in polar coordi- 
nates, which quickly destroys the gas disk during the nu- 
merical simulations. The exponent a(r) decreases linearly 
with radius r. More specifically, we choose a(r = kpc) = 
2.3, and a(r = 17 kpc) = 0.6. This specific form of C(r) 
allows us to produce spiral gravitational potentials with dif- 
ferent shapes and maximum amplitudes. The resulted ratio 
f3(r) of the maximum non-axisymmetric gravitational accel- 
eration due to spiral arms [(d$ sp /dr) 2 + (r~ 1 <9$ sp / 1 d4>) 2 ] 1 / 2 
(at a given radial distance r) versus the gravitational ac- 
celeration due to both the dark matter halo d^^/dr and 
axisymmctric stellar disk <9$ st /<9r is shown in Fig. Q] by the 
solid line. We note that ji never exceeds 19%. 

The measurements of iKuzio de Narav et al.l (|2006f) in- 
dicate that the gas velocity dispersions of LSB galaxies lie 
in the range between 6 km s _1 and 10 km s _1 and are con- 
sistent with the typical dispersions for the gas component 
of HSB galaxies. Therefore, we assume that the gas disk 
of our model galaxy is initially isothermal at a tempera- 
ture T = 10 4 K. LSB galaxies have same atomic hydro- 
gen content as their HSB counterparts bu t generally lower 
surface densit i es an d extended gas disks (jPickerind . Il999t 
Ide Blok et all |1996| ) . We assume that the gas disk has an 
exponentially declining density profile 




6 8 10 12 

Radial distance (kpc) 

Fig. 1. (Solid line.) Ratio (3(r) of the maximum non- 
axisymmetric gravitational acceleration due to spiral arms 
[(d^ sp /dr) 2 + (r^ 1 d$ sp / d^) 2 ] 1 / 2 vesus gravitational ac- 
celeration due to both the dark matter halo <9$h/<9r and 
axisymmctric stellar disk 3$ st /9r. (Dashed line.) Vertical 
scale height Z of our model gas disk as a function of radius 
in the beginning of numerical simulations (t = 0). 
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Fig. 2. The initial rotation curve (solid line) and gas surface 
density distribution (dashed line). 



with the central surface density S g o = 6.5 M pc -2 and 
radial scale length r g = 30 kpc. The gas disk extends 
to 17 kpc in our model, but star formation is confined 
to the inner 15 kpc. The outermost 2 kpc are kept for 
numerical reasons. Gas disks of such extent (~ 1 5 kpc ) 
are typical for LSB galaxies (jvan den Hoek etaTl |2000| ). 
The mean surface density of our model disk is (S g ) = 
4.9 M Q pc -2 , which is comparable to a median HI surface 
de nsity Shi = 3.8 Mp, pc~ 2 of a sample of LSB galaxies 
in Ivan den Hoek et"~atl (|2000f ) . The mass of our model gas 
disk in the inner 15 kpc is M g = 3.3 x 10 9 M Q . Hence, the 
gas disk contains only a small fraction (~ 15%) of the total 
mass in the computational domain and, to a first approxi- 
mation, we can neglect its self-gravity. 

The initial rotation curve of the gas disk is calculated by 
solving the steady-state momentum equation for the tan- 
gential component of the gas velocity 



£ g = E g0 exp(-r/V g ), 



(6) 



RT (9£ s <9$ st d$ h 



Sg/i dr 



dr 



dr 



(7) 
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Table 1. Main structural properties of our model galaxy 





mass 


scale length 


central density 


Stellar disk 


2.6 x 10 a 


4 


30 


gas disk 


3.3 x 10 9 


30 


6.5 


halo 


2.0 x 10 10 


5.7 


6.0 x 10~ 3 



All masses are in Mq, scale lengths in kpc, and densities in 
Mq pc~ 2 (gas and stellar disks) and Mq pc~ 3 (halo). All 
masses are calculated inside 15 kpc radius. 



where /i = 1.2 is the mean molecular weight for a gas of 
atomic hydrogen (76% by mass) and atomic helium (24% 
by mass) and R is the universal gas constant. The initial 
radial gas density profile and rotation curve are shown in 
Fig. [U by the dashed and solid lines, respectively. It is evi- 
dent that our model LSB galaxy is characterized by a rising 
rotation curve with a maximum rotation velocity of approx- 
imately 80 km s -1 , in agreement with the measured shapes 
of the rotation curve in UGC 5750 and other LSB galaxies 
iKuzio de Narav et all (|2006f) . The main properties of our 
model galaxy are summarized in Table [TJ 

3.2. Basic equations 

We use the thin-disk approximation to compute the long- 
term evolution of our model LSB galaxy. It is only the 
gas disk that is dynamically active, while both the stel- 
lar disk and dark matter halo serve as a source of exter- 
nal gravitational potential. In the thin-disk approximation, 
the vertical motions in the gas disk are neglected and the 
basic equations of hydrodynamics are integrated in the z- 
direction from z = —Z(r, <fi, t) to z = +Z(r, 4>, t) to yield 
the following equations in polar coordinates (r, <j)) 

dt 



V p • («pSg) — —0RM SsFR, 



(8) 



ds p 
dt 



de 
dt 



+ V p • (s p ■ v p ) = -V p p - E g Vp$ h 
- S g Vp($ sp + $st) - v p p R M S SF r, 

v p • (v p e) = -p(v ■ v) + r sn + 2Z(r cr + r bg ) 



(9) 



- 2ZA — v'pPrm Ssfr- (10) 



In the above equations, v p = fv r + 4>v^, is the gas velocity 
in the disk plane, s p = S g i>p is the momentum per surface 
area in the disk plane, V p = fd/dr + (j)r~ 1 d/d<t) is the 
gradient along the planar coordinates of the disk, p is the 
vertically integrated gas pressure, e is the internal energy 
per surface area, and Z is the vertical scale height of the gas 
disk. The basic equations of hydrodynamics are modified to 
include the effect of star formation, where Esfr is the star 
formation rate (SFR) per unit area and /3rm = 0.42 is the 
fraction of stellar mass tha t ends up as an inert remnant 
(jKoppen fe Arimotol . I1991D . Equations ([g])-([T0j) are closed 
with the equation of state of a perfect gas p = (7 — l)e and 
7 = 5/3. 

The volume cooling rate A = A (erg cm 3 s as a 
function of g as temperature T i s calc ulated from the cool- 
ing curves of IWada fc Norrrianl (|200ll ) and is linearly inter- 
polated for metallicities between 10~ 4 and 1.0 of the solar. 



The heating of gas is provided by supernova energy release 
r sn , cosmic rays F cr , and background ultraviolet radiation 
field Tt,g . The supernova energy release (per unit time and 
unit surface area) is defined as 



r — 

1 sn — 



e sn SsFR/8MQ m ^ m ) rfm 

/™ up m^(m)dm' 



(11) 



e sn = 10 51 erg is the energy release by a single supernova, 
and m sn = 15 Mq is the mean mass of massive stars capa- 
ble of producing supe rnova explosio ns. We use the Kroupa 
initial mass function (jKroupal l200ll ) defined as 



£(m) 



4m" 1 ' 3 if to < 0.5 M Q 
BtoT 2 - 3 ifTO>0.5M P 



1 



(12) 



with the lower and upper cutoff masses m\ ow = 0.1 Mq and 
TO up = 100 Mq, respectively. The ratio of the normalization 
constants A/B = 2 is found using the continuity condition 
at m = 0.5 Mq. We adopt the cosmic ray heating per unit 
time 



T cr = 10~ 27 n g (erg cm" 3 s" 1 ) , 



(13) 



where n g = £ g / (2Z fxmn) is the volume number density 
of gas. Heating due to the background ultraviolet radia- 
tion field is supposed to balance cooling and help sustain 
a gas temperature T = 10 4 K in the unperturbed by star 
formation gas disk of our model galaxy. Therefore, Tb g is 
determined in the beginning of numerical simulations from 
the following relation Fb g ?i g = A, where A and n g are those 
of the initial gas disk configuration, and is kept fixed after- 
wards. 

The volume heating and cooling rates in equation (jlpp 
are assumed to be independent of z-direction and are mul- 
tiplied by the disk thickness (2Z) to yield the correspond- 
ing vertically integrated rates. The vertical scale height 
Z = Z(r,(j),t) is calculated via the equation of local pres- 
sure balance in the gravitational field of the dark matter 
halo, stellar, and gas disks 



= 2 



Ps [9% 



^halo 



)dz, 



(14) 



where p g c 2 is the gas pressure in the disk midplane, gf as , 
gf, and are the vertical gravitational accelerations 

due to self-gravity of a gas layer, stellar layer, and grav- 
itational pull of the dark matter halo, respectively, and 
c s = 7(7 ~ l) e /^g i s tne square of the sound speed. 
The details of this method are given in the Appendix. 
We plot Z as a function of radius by the dashed line in 
Fig. [T] for our model gas disk in the beginning of numeri- 
cal simulations. It is evident that the disk is rather thick, 
which is indeed expected for an intermediate-mass galaxy 
(~ a few x 10 10 Mq) considered in the present paper. In 
practice, we do not let the vertical scale height to drop be- 
low 100 pc, since such low values of Z are not expected in 
LSB galaxies and may lead to overcooling of our model gas 
disk. It should be noticed that in our numerical simulations 
the energy release by supernova explosions is localized in 
the disk due to the adopted thin-disk approximation. In 
reality, some portion of the supernova energy may be lost 
to the intergalactic medium if the vertical scale height of a 
gas disk is sufficiently small and the energy release by mul- 
tiple supernova explosions is sufficiently large. Given that 
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intermediate-mass LSB galaxies are characterized by thick 
gas disks and low rates of star formation (~ 0.1 Mq yr -1 ), 
we do not expect this blowout effect to considerably influ- 
ence the dynamics of our model gas disk. 

3.3. Instantaneous recycling approximation 

The temporal and spatial evolution of oxygen in our model 
galaxy is computed adopting the so-called instantaneous 
recycling approximation, in which an instantaneous enrich- 
ment by freshly synthesized elements is assumed. This ap- 
proximation is valid for the oxygen enrichment, because 
most of oxygen is produced b y short-lived, massive stars . 
We adopt the oxygen yields of lWooslev fc Weaver! (|1995f ). 
We assume that oxygen is collisionally coupled to the gas, 
which eliminates the need to solve an additional momen- 
tum equation for oxygen. The continuity equation for the 
oxygen surface density S ox modified for oxygen production 
via star formation reads as 



dt 



(15) 



The oxygen production rate per unit area by supernovae is 
calculated as 



/i2aL P( m ) £(m)dm 



J SFR" 



/12M, 

T 

■J n 



(16) 



where p( m) is the mass of oxygen released by a star with 
mass m ' (fWooslev fc Weaveii Tl99 5h. We assume that the 
initial oxygen abundance of our model galaxy is [O/H] = 
log 10 (£ ox /£ g ) - log 10 (E ox /Sg) Q = -4, where the ratio 
(E ox /S g ) Q in the Solar neighbourhood is set equal to 
7.56 x 1CT 3 . 

3.4. Code description 

An Eulerian finite-difference code is used to solve equa- 
tions (fg|)- (fTU)) and (fT5|) in polar coordinates (r, (f>). The ba- 
sic algorithm of the code is similar to that of the commonly 
used ZEUS code. The operator splitting is utilized to ad- 
vance in time the dependent variables in two coordinate 
directions. The advection is treated using the consistent 
transport method and the van Leer interpolation scheme. 
There are a few modifications to the usual ZEUS method- 
ology, which are necessitated by the presence of star for- 
mation and cooling and heating processes in our model hy- 
drodynamic equations. The total hydrodynamic time step 
is calculated as a harmonic average of the usual time step 
due to the Courant-Fricdrichs-Lewy condition and the star 
formation time step defined as isF = O.lEg/SgpR. The in- 
ternal energy update due to cooling and heating in equa- 
tion (|10p is done implicitly via Newton-Raphson iterations, 
supplemented by a bisection algorithm for occasional zones 
where the Newton-Raphson method does not converge. 
This ensures a considerable economy in the CPU time. A 
typical run takes about four weeks on four 2.2 GHz Opteron 
processors. The resolution is 500 x 500 grid zones, which are 
equidistantly spaced in both coordinate directions. The ra- 
dial size of the grid cell is 34 pc and the radial extent of 
the computational area is 17 kpc. However, star formation 
is confined to only the inner 15 kpc. The code performs 
well on the angular momentum conservation problem (sec 



e.g. lVorobvov fc Basul . l2006l ). This test problem is essential 
for the adequate modelling of rotating systems with radial 
mass transport. 

3.5. Sporadic star formation 

Low gas surface densities (|de Blok et all . 119961 : iPickerinel 
Il999f ) , which are usually below the critical value determined 
by the Toomre criter i on, an d low total star formation rates 
van den Hoek etafl |200(I) argue against large-scale star 
formation as a viable scenario for LSB galaxies. Despite 
the low gas surface densities, there is evidence for a no- 
ticeable young stellar population as i ndica ted by the blue 
colors of most LSBs (|de Blok et all [19951). Ac c ording to 
Ide Blok et ail lj!995f) and iGerretsen fc de Blokl l)1999h . a 
sporadic star formation scenario (i.e. small surges in the 
SFR, either superimposed on a very low constant SFR or 
not) can best explain the observed blue colors and evidence 
for young stars. 

There is little known about physical processes that con- 
trol sporadic star formation in LSB galaxies. The recent 
uncalibrate d Ha imaging of a sample of LSB galaxies by 
lAuld et ail |2006:) and the observational data presented 
here (see Sect.HJ) suggest that star formation is highly clus- 
tered. Current star formation is localized to a handful of 
compact regions. There is little or no diffuse Ha emission 
coming from the rest of the galactic disk. Motivated by this 
observational evidence, we adopt the sporadic scenario for 
star formation in our model LSB galaxy. 

We use a simple model of sporadic star formation in 
which the individual star formation sites (SFSs) are dis- 
tributed randomly throughout the galactic disk and the star 
formation rate per unit area in each SFS is determined by 
a Schmidt law 



J SFR 



= «SF 



M Q pc- 



1.5 



if T < T CA 



(17) 



M Q yr _1 kpc 

where q;sf = 6 x 10~ 4 is a proportionality constant and 
T cr = 2 x 10 4 K. If the gas temperature in a randomly 
chosen SFS exceeds the critical value, the SFS is rejected. 
The adopted value of T cr = 2 x 10 4 K differs from typical 
temperatures of star-forming molecular clouds: T ^ 30 — 
50 K. They are however more appropriate for identification 
of star-forming regions on ~ 34 pc scales that are resolved 
in the present numerical simulations. A Schmidt law with 
index n ~ 1.5 would be expected for self-gravitating disks, 
if the SFR is equal to the ratio of the local gas volume 
density (p) to the free-fall time (oc p~ - 5 ), all multiplied 
by some star formation efficiency e. The efficiency e is a 
measure of the fraction of the gas mass converted into stars 
before the clouds are disrup t ed. It can be shown that e is 
related to cxsf as (IVorobvovL[2003h 



q sf = 0.12 e 



pc 



-0.5 



(18) 



This means that the adopted star formation efficiency in 
our model ranges between 0.045 and 0.15 for possible values 
of the vertical scale height Z = 100 — 1500 pc. 

The number of star formation sites at any given time 
can be roughly determin ed using the uncal ibrated Ha im- 
ages of LSB galaxies by (|Auld et all |2006[ ). These images 
highlight only a handful of small regions that have re- 
cently formed stars. There is little diffuse Ha emission, 
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contrary to the emission in K-. R-, and B-bands which (as 
a rule) show a noticeable diffuse component. The number of 
HII regions within an indi vidual LSB galaxy identified by 
van den Hoek et~aTl (|2000h ranges between 4 and 26, with 
a median number of HII regions equal to 11. Therefore, the 
number of SFS, iVgps, in our numerical simulations is set 
to approximately 15. If we further assume that both N^fs 
and a typical area occupied by an individual star forma- 
tion site Ssfs do not change significantly with time, this 
would yield a near-constant mean SFR over the lifetime 
of our model galaxy. However, there is evidence that star 
formation in blue LSB galaxies does not proceed at a near- 
co nstant rate. Our own nu merical simulations and modeling 
bv lZackrisson et aTl (|2005f ) indicate that the SFR should be 
declining with time to reproduce the observed Ha equiv- 
alent widths in LSB galaxies. Therefore, we assume that 
the typical area occupied by a SFS is a decreasing func- 
tion of time of the form Ssfs = Ssfso ex P(~ V t sfr), where 
Ssfso = 450x450pc is the area at t = and tsfr = 1.3 Gyr 
is a characteristic time of exponential decay in the star for- 
mation rate. 



The duration of star formation in an individual star 
formation site (tsfs) is poorly known. We use the theoret- 
ical and empirical knowledge gained by analyzing similar 
processes in the Solar neighbourhood. Observations of the 
current star formation activity in nearby molecular cloud 
complexes indicate that t he ages of T Tauri sta rs usually 
fall in the 1-10 Myr range (jPalla fe Stahlerl . l2002D . Clusters 
that are 10 Myr old often have little associated gas remain- 
ing, implying that the process shuts off by then. However, 
the volume gas density in disks of LSB galaxies is expected 
to be lower (on average) than that of the Solar neighbour- 
hood, which may prolong the process of gas conversion into 
stars. It should also be noted that individual star forma- 
tion sites in LSB galaxies have characteristic sizes that are 
much larger than those of an individual giant molecular 
cloud. Otherwise, we would have to assume unrealistically 
high star formation efficiency in LSB galaxies to retrieve the 
observed star formation rates of order 0.1 Mq yr _1 . This 
implies a presence of local triggering mechanisms within an 
individual SFS that should also act to prolong the star for- 
mation activity. From a theoretical background one can ar- 
gue that the duration of star formation should be limited by 
a feedback action from newly born massive stars and super- 
nova explosions. The lifetime of a least massive star capable 
of producing a supernova is approximately 40 Myr. Hence, 
the actual lifetime of an individual SFS should lie in be- 
tween 10 Myr and 40 Myr. Using spectro- photometric evo- 
lution models, Ivan den Hoek et al.l (|2000l ) have also found 
that tsfs = 10 — 50 Myr reproduce well the observed (R—I) 
colors and /-band luminosities of LSB galaxies. In our nu- 
merical simulations we assume that tsfs = 20 Myr. Thus, 
each twenty million years we initialize a new set of SFSs 
(•^sfs ~ 15) and randomly distribute them in the inner 
15 kpc of the galactic disk. Each SFS is assigned a velocity 
equal to that of the local gas velocity upon its creation and 
is evolved collisionlessly in the total gravitational potential 
of the stellar disk and the dark matter halo using the fifth- 
order Rungc-Kutta Cash-Karp integration scheme. The old 
SFSs are rendered inactive and their further evolution is not 
computed. 
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Fig. 3. Reprocessed and continuum subtracted images of 
the LSB galaxies UGC 9024 (upper left), UGC 12845 (up- 
per right), UGC 5675 (lower left), and UGC 6151 (lower 
right). 



4. Observations 

There is still just a small number of (nearby) LSB galaxies 
with published Ha images, which would allow the study of 
the distribution and properties of their star-formation re- 
gions. Ha images of about 1 8 LSB galaxies are avai lable in 
McGaugh fc Bothunl dl99l; iMcGaugh et all (jl995h . while 
de Blok fc van der Hulstl ( 19981 ) only published the spec- 
tra and some broad band images. A new sa mple of 6 LSB 
galaxies was presented in lAuld et ail (|2006[ ) . 

From the image data set (McGaugh, priv. comm.) we 
selected 4 LSB spiral galaxies based on their B band mor- 
phology. We selected the galaxies to be close to face-on and 
showing different disk surface brightness and spiral mor- 
phology. This resulted in UGC 9024, UGC 12845, UGC 
5675, and UGC 6151. A different version of the Ha images 
of UGC 902 4, UGC 5675, and UG C 6151 was previously 
published in lMcGaugh et all (|l995l ). 

For the presentation in Figure[3]we improved on the cos- 
mic rays rejection and correction for detection blemishes of 
the already flat-fielded images. The generation of the con- 
tinuum correcte d Ha images fo llowed closely the schema 
described e.g. in iBomansl (|l997h . We realigned the Ha and 
continuum images carefully using the foreground stars visi- 
ble in both images to sub-pixel accuracy and subtracted the 
scaled continuum image. We tested different scaling factors 
and defined an optimal subtraction based on the disappear- 
ance of the stellar disk of the LSB galaxy. These factors also 
gave the smallest foreground star residuals and is consistent 
with the scaling factor calculated by measuring the flux of 
several foreground stars in both images. 

The distribution of the HII regions in all 4 galaxies in 
Figure[3]is consistent with the predicted pattern of sporadic 
star formati on. The Ha i mage s of 6 LSB galaxies recently 
published by lAuld et alj (|2006l ) sh ow the same behavior, a s 
do the other galaxies published in iMcGaugh et al.l (|l995h . 
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5. Numerical results 

5.1. Gas density and temperature distribution 

A model galaxy described in Section [3] will be referred here- 
after as model 1. In a due course, we will also introduce 
other models that will have identical structural properties 
such as those summarized in Tabic [TJ but different IMFs and 
star formation rates. We start our numerical simulations by 
setting the equilibrium gas disk and slowly introducing the 
non-axisymmetric gravitational potential of the stellar spi- 
ral arms. Specifically <& S p is multiplied by a function 
which has a value of at t = and linearly grows to its max- 
imum value of 1.0 at t > 200 Myr. It takes another hundred 
million years for the gas disk to adjust to the spiral gravi- 
tational potential and develop a spiral structure. Figure [4] 
shows the gas surface density distribution (in log units) 
in the 1-Gyr-old (upper left), 2-Gyr-old (upper right), 5- 
Gyr-old, (lower left), and 12-Gyr-old (lower right) disks. It 
is seen that the gas disk develops a strong phase separa- 
tion with age. While the initial gas disk has surface densi- 
ties between 3.5 M Q pc -2 and 6.5 M pc -2 , the 1-Gyr-old 
disk features some dense clumps with surface densities in 
excess of 100 M© pc~ 2 and temperatures of order 100- 
200 K. The number of such clumps increases with time and 
they are getting stretched into extended filaments and arcs 
by galactic differential rotation. Observational detection of 
such features could be a powerful test on the sporadic na- 
ture of star formation in LSB galaxies. The gas response to 
the underlying spiral stellar density wave is rather weak in 
the early disk evolution and can hardly be noticed in the 5- 
and 13-Gyr-old disks. It implies that mixing due to spiral 
stellar density waves is inefficient in the old disk. 

The phase separation of the gas disk can be nicely il- 
lustrated by calculating the gas mass per unit tempera- 
ture band at different evolutionary times. Figure [3] shows 
the gas mass per Kelvin as a function of temperature at 
t = 2 Gyr (solid line), t = 5 Gyr (dashed line) and 
t = 13 Gyr (dash-dotted line). There are two prominent 
peaks that manifest the separation of the gas into two 
phases - a cold gas phase with peak temperature about 
T = 50 — 100 K and a warm gas phase with peak tem- 
perature about T = (1.0 - 1.5) x 10 4 K. The positions of 
the peaks drift apart as the gas disk evolves, indicating 
that the cold phase becomes colder and the warm phase 
becomes warmer with time. However, the gas temperature 
never drops below 35 K. It is also evident from the ampli- 
tudes of the peaks that the cold phase accumulates gas mass 
with time, which is explained by increased cooling due to 
ongoing heavy element production. A similar phase separa- 
tion was also obtained in nu merical simulations of LSBs by 
iGerretsen fc de Blokl (| 19991) . The fact that there is little gas 
with temperature below 35 K implies that most of the cold 
phase may be in the form of atomic rather than molecular 
hydrogen, but more accurate numerical simulations with a 
chemical reaction network and higher numerical resolution 
are needed to prove this hypothesis. 

Figure [5] indicates that model 1 has no hot gas phase 
with temperatures of order 10 6 K, which is a consequence 
of a relatively low SFR and moderate numerical resolution. 
The energy injection by star formation of such a low rate 
is not enough to convert all gas in a computational cell 
(~ 34 x 34 pc) to a hot phase. To remind the reader, we 
solve the equations of hydrodynamics for a single fluid. As a 
consequence, temperature of all gas in a computational cell 
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Fig. 4. Gas surface density distribution in the 1-Gyr-old 
(upper left), 2-Gyr-old (upper right), 5-Gyr-old (lower left), 
and 12-Gyr-old (lower right) disks. Dark regions are the 
densest and coldest. The scale bar is Mq pc -2 (log units). 
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Fig. 5. Gas mass per Kelvin versus temperature in disks 
of distinct age as indicated in the legend. The two peaks 
manifest the phase separation. Note that there is little gas 
with temperatures below 40 K. 



is characterized by a specific value, there is no gas phase 
separation within an individual cell. This, however, does 
not mean that the hot gas is not present at all in LSB 
galaxies, we simply may need a much higher numerical res- 
olution to resolve it. We do start seeing the hot phase if we 
increase the SFR in our mode l by a factor of several. We 
note that IGerretsen fc de B lok (1999]) have also found little 
hot gas in their numerical simulations. 
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Fig. 6. Star formation rate versus time in model 1. The 
solid and thick dashed lines show the SFR averaged over 
20 Myr and 1.0 Gyr, respectively. Note large amplitude 
fluctuations in the 20-Myr-averaged SFRs around the 1- 
Gyr-averaged values. 
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Fig. 7. Radial profiles of [O/H] obtained by making a nar- 
row cut through the galactic disk at an azimuthal angle 
= 0°. The age of the disk is indicated in each frame. 



Figure [S] shows the SFR as a function of time in model 1. 
In particular, the solid line presents the SFR averaged over 
the past 20 Myr, these values are later used to construct 
Ha equivalent widths. In addition, the thick dashed line 
plots the SFR averaged over a 1 Gyr interval centered on 
the current time (running average method). It is seen that 
the 20-Myr-averagcd SFRs have considerable fluctuations 
around the 1-Gyr-averaged values. These fluctuations have 
a characteristic time scale of a few hundreds of Myr and 
amplitudes that may exceed 0.1 Mq yr . The existence 
of such bursts of star formation is an essential ingredient 
for successful modelling of blue colors of at least some LSB 
galaxies (e.g. Ivan den Hoek et ail . l2000Ti . Similar fluctua- 
tions in the SFR were also observed in numerical simula- 
tions bv iGerretsen fc de Blokl (|1999t ). The 1-Gyr -averaged 
SFR quickly grows to a maximum value 0.26 Mq yr -1 at 
t = 2.6 Gyr and subsequently declines to 0.08 Mq yr -1 
at t = 13.0 Gyr. Typical star formation rates in LSB 
galaxies range from 0. 02 to 0.2 Mq yr" 1 , with a median 
SFR of 0.08 M Q yr" 1 (jvan den Hoek et al.l . [200l . Model 1 
yields SFRs that are consistent with the inferred present- 
day rates. 

5.2. Radial variations in the oxygen abundance 

We calculate the oxygen abundance [O/H] = 
log 10 (E ox /E g ) — log 10 (7.56 x 10~ 3 ) in our disk from 
the model's known surface densities of gas (S g ) and oxy- 
gen (S ox ). We find that the spatial distribution of oxygen 
in the 1-Gyr-old gas disk has a pronounced patchy ap- 
pearance - local regions with enhanced oxygen abundance 
adjoin local regions with low oxygen abundance. The local 
contrast in [O/H] may be as large as 1.0 dex or even more. 
As the galaxy evolves, the spatial distribution of oxygen 
becomes noticeably smoother. The smearing is caused by 
stirring due to differential rotation, radial migration of gas 
triggered by the non-axisymmetric gravitational potential 
of stellar spiral arms, and energy release of supernova 
explosions. 



To better illustrate the spatial variations in the oxygen 
distribution in model 1, we show in Fig.[7]thc radial profiles 
of the oxygen abundance obtained by making a narrow (~ 
50 pc) radial cut through the galactic disk at = 0°. The 
numbers in each frame indicate the age of our model galaxy. 
The most prominent feature in Figurc[7]is a sharp decline in 
a typical amplitude of the radial fluctuations in [O /H] . For 
instance, the radial distribution of the oxygen abundance 
in the 1-Gyr-old disk is very irregular and shows radial 
fluctuations with amplitudes as large as 1.0 dex and even 
more. The fluctuations have a short characteristic radial 
scale of order 1 — 2 kpc. The fluctuation amplitudes have 
noticeably decreased in the 2-Gyr-old disk, which now has 
a maximum amplitude of order 0.5 dex. The 5-Gyr-old disk 
(and older disks) feature only low-amplitude fluctuations 
in the oxygen abundance with amplitudes rarely exceeding 
0.2 dex. 

Figure [7] reveals the absence of any significant radial 
abundance gradients in our mo del disk. This appears to 
be ty pical for LSB galaxies (e.g. Ide Blok fc van der Hulstj . 
Il998l ) . The lack of radial abundance gradients is most likely 
caused by a spatially and temporally sporadic nature of 
star formation, i.e., LSB galaxies must have evolved at the 
same rate over their entire disk. Well-known star formation 
threshold criteria, such as those based on the gas surface 
density threshold (Toomre c riterion) or the r ate of shear 
(|Martin fc Kennicuttj (|200lh ; IVorobvovi (|2003l )L would in- 
crease the likelihood for star formation in regions of higher 
density and lower angular velocity, which are usually the in- 
ner galactic regions. As a result, negative radial abundance 
gradients may develop with time. The fact that we do not 
see radial abundance gradients in LSB galaxies implies that 
star formation has no threshold barrier. 

To analyze the properties of the oxygen abundance fluc- 
tuations in disks of different age, we calculate the fluctua- 
tion spectrum defined as the (normalized) number of fluctu- 
ations versus the fluctuation amplitude A. The latter is cal- 
culated in the following manner. We scan the disk in the ra- 
dial direction along each azimuthal angle ((f>) of our numer- 
ical polar grid for the adjacent minima and maxima in the 
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Fig. 8. Oxygen abundance fluctuation spectrum defined as 
the normalized number of radial oxygen abundance fluctu- 
ations versus fluctuation amplitude (in dex). Lines of dis- 
tinct color corresponds to different numerical models and 
lines of distinct type corresponds to different galactic ages 
as indicated. See the text for more details. 

oxygen abundance distribution. The absolute value of the 
difference between the adjacent minimum and maximum 
(or vice versa) gives us the radial amplitude of a fluctua- 
tion on spatial scales of our numerical resolution, ~ 40 pc. 
The calculated fluctuation amplitudes are distributed in 40 
logarithmically spaced bins between 0.1 dex and 3.0 dex 
and then normalized by the total number of fluctuations in 
all bins. The resultant fluctuation spectrum has the mean- 
ing of the probability function ^(A) — the product of the 
total number of fluctuations and F(A) yields the number 
of fluctuation with amplitude A. 

Figure [5] shows the probability function J- (A) derived 
for four models detailed below. Hereafter, different line 
types in Fig. [5] will help to distinguish between galaxies 
of distinct age and different line colors will correspond to 
different numerical models. Black lines present J- (A) in 
model 1 for the 1-Gyr-old disk (solid line), 2-Gyr-old disk 
(dashed line), 5-Gyr-old disk (dash-dotted line), and 12- 
Gyr-old disk (dash-dot-dotted line). It is evident that the 
probability functions for disks of different age intermix for 
T{A) ^ 0.1 but appear to be quite distinct for T(A) < 0.05. 
Can this peculiar feature of J- (A) be used to constrain the 
age of LSB galaxies? The extent to which oxygen is spa- 
tially mixed in the galactic disk may depend on the model 
parameters such as the strength of the spiral density wave, 
shape of the IMF, and details of the energy release by SNc. 
In order to estimate the dependence of J- {A) on these fac- 
tors, we run two test models. In particular, model Tl has 
no spiral density waves and model T2 is characterized by 
a factor of 2 lower energy input by SNe (latter 'T' stands 
for 'test' to distinguish from regular models). Model T2 at- 
tempts to mimic the case when part of the SN energy is lost 
to the extragalactic medium via SN bubbles. It is also im- 
portant to see the effect that a different IMF may have on 
the probability function. Therefore, we have also introduced 
model 2, which is characterized by the Salpeter IMF with 
the upper and lower cut-off masses of 0.1 M Q and 100 M Q . 
Model 2 has the same structural parameters as model 1 
except for tsfr which is set to 18.5 Gyr. The parameters 



of models 1 and 2 are detailed in Table [U The probabil- 
ity functions corresponding to model Tl, model T2, and 
model 2 arc shown in the top panel of Fig. [8] by the red, 
blue, and cyan lines, respectively. The meaning of the line 
types (i.e., solid, dashed, etc.) is indicated in the figure. 

A noticeable scatter in the probability functions of dis- 
tinct age is seen among different models. In particular, the 
probability functions of galaxies of 5-12 Gyr age intermin- 
gle and thus cannot be used for making age predictions. 
However, it is still possible to differentiate between (1-2)- 
Gyr-old and (5-12)-Gyr-old galaxies based on the slope of 
J- (A) . Young galaxies of 1-2 Gyr age have fluctuation spec- 
tra with a much shallower slope than (5-13)-Gyr-old galax- 
ies, meaning that the former feature radial oxygen abun- 
dance fluctuations of much larger amplitude than the lat- 
ter. A similar behaviour of J~(A) is found when we decrease 
the spatial resolution of the fluctuation amplitude measure- 
ments from 40 pc to several hundred parsecs (more appro- 
priate for modern observational techniques). We conclude 
that measurements of the radial oxygen abundance fluc- 
tuations may be used only to set order- of -magnitude con- 
straints on the age of blue LSB galaxies. 

We are aware of only one successful measurement of the 
radial abundance gradients in t hree LSB galaxies done by 
Ide Blok fc van der Hulstl (|T998). Unfortunately, the num- 
ber of resolved HII regions and oxygen abundance measure- 
ments in each galaxy is too low (< 10) to make a mean- 
ingful comparison with our predictions. With such a small 
sample, the probability function would be confined in the 
0.1 < J- {A) < 1.0 range, which is not helpful for making 
age predictions due to degeneracy of J- (A) in this range. 

6. Constraining the age of blue LSB galaxies 

In this section, we investigate if the time evolution of the 
mean oxygen abundance, complemented by Ha emission 
equivalent widths and optical colors, can be used to better 
constrain the age of blue LSB galaxies. We consider several 
new models (in addition to our models 1 and 2) that have 
different rates of star formation. 

6.1. Time evolution of the oxygen abundance 

There have been several efforts in the past to measure 
the oxygen abundances in HII regions of blue LSB galax- 
ies. For ins tance, oxygen abundances in a sample of 18 
galaxies in iKuzio de Narrav et alj (|2004f ) have a median 
value ([O/H]) = —0.95. Their smallest and largest mea- 
sured abundances (using the equivalent width method) are 
JO /HI = -12 (F611-1) a nd [O /H] = -0.04 (UGC 5709). 
Ide Blok fc van der Hulstl (|1998f ) have observed 64 HII re- 
gions in 12 LSB galaxies and derived oxygen abundances 
that peak around [O/H] = —0.80. The lowest abundance in 
their sample is [O/H] = —1.1. These data demonstrate that 
the present day mean oxygen abundance in LSB galaxies is 
about {[0/H])lsb = —0.85 dex, albeit with a wide scatter. 
We plot the minimum and maximum oxygen abundances 
found in blue LSB galaxies by horizontal dotted lines in 
Fig. [9l while the mean oxygen abundance for LSB galaxies 
is shown by the vertical dash-dot-dotted line. 

The amount of produced oxygen and the oxygen abun- 
dance are expected to depend on the adopted rate of star 
formation and on the details of the initial mass function. 
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Table 2. Summary of galactic models 



model 


IMF type 


upper/lower cutoffs 


ISFR 


Q?SF 




min/max SFR 


1 


Kroupa 


0.1 - 100 M© 


13.0 


6 x 10" 


4 


0.075/0.26 


1L 


Kroupa 


0.1 - 100 Mq 


13.0 


4 x 10" 


-4 


0.085/0.17 


2 


Salpeter 


0.1 - 100 Mq 


18.5 


6 x 10" 


-4 


0.045/0.20 


2L 


Salpeter 


0.1 - 100 Mq 


18.5 


4 x 10" 


-4 


0.050/0.14 



Star formation rate is in Mq yr and characteristic time of SFR decay tsfr is in Gyr. 
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Fig. 9. Mean oxygen abundance as a function of galac- 
tic age. Each line type corresponds to a particular nu- 
merical model as indicated in the figure. The horizon- 
tal dotted lines delineate t he range of oxygen abundances 
found in blue LSB galaxies (Ide Blok fc van der Hulstl . Il998t 
iKuzio de Narrav et all I2004 ) and the horizontal dashed- 
dot-dotted line shows the mean oxygen abundance derived 
from a large sample of LSB galaxies. See the text and 
Table [3 for more details. 



Most of oxygen is produced by intermediate to massive stars 
and both the slope and the upper mass cutoff are those pa- 
rameters of the IMF that may influence the rate of oxygen 
production. In this paper, we study IMFs of different slope, 
leaving cut-off IMFs for a subsequent study. 

The mean oxygen abundances ([O/H]) obtained in 
model 1 are presented in Fig. \§\ by the thick solid line. The 
mean values are calculated by summing up local oxygen 
abundances in every computational cell between 1 kpc and 
14 kpc and then averaging over the corresponding number 
of cells. The dashed line shows ([O/H]) as a function of 
galactic age in model 2. In addition, we have introduced 
two more models that have structural parameters identical 
to those of models 1 and 2 (see Table [T]) but are character- 
ized by lower SFRs. In particular, models 1L and 2L (thin 
solid line/dash-dotted line) feature on average a factor of 
1.4 lower rates of star formation than models 1 and 2. The 
parameters of models 1L and 2L, as well as the minimum 
and maximum SFRs, are listed in Tabled It is evident that 
models with Salpeter IMF yield lower ([O/H]) than the cor- 
responding models with Kroupa IMF. With the lower and 
upper mass limits being identical, the Sapleter IMF may be 
regarded bottom-heavy in comparison to that of Kroupa. 
A bottom-heavy IMF has a lower relative number of mas- 
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Fig. 10. Ha equivalent widths (left column) and global 
B — V colors (right column) for model 1 (top row) and 
model 2 (bottom row). The horizontal dotted lines delineate 
a range of observed val ues i n blue LSB galaxies acc ording 
to Ide Blok et all (| 19951) and IZackrisson etall (|2005l ). The 
thick dashed lines show EWKa averaged over 1 Gyr. The 
arrow shows a (possible) reddening due to dust extinction 
for a fiduc i al val ue of the color excess E(B — V) = 0.1 
riMcGaueE 11991 . 



sive stars and this fact accounts for a lower rate of oxygen 
production in models with the Salpeter IMF. Furthermore, 
models with lower SFRs yield lower mean oxygen abun- 
dances. 

There are several interesting conclusions that can be 
drawn by analyzing Fig.[9J First, our models can reproduce 
the entire range of observed oxygen abundances found in 
LSB galaxies (delimited by horizontal dotted lines). Second, 
for a given value of ([O/H]), the predicted ages show a large 
scatter (that increases along the sequence of increasing oxy- 
gen abundances), making mean oxygen abundances alone a 
poor age indicator, especially for metal-rich LSBs. Finally, 
the minimum oxygen abundance found among LSB galax- 
ies (bottom dotted line in Fig [9]) can be attained after just 
1-3.5 Gyr of evolution, suggesting a possibility of rather 
young age for LSB galaxies. We will return to this issue in 
some more detail in the following section. 

6.2. Ha equivalent widths and B-V integrated colors 

In this section we use population synthesis modelling to cal- 
culate the integrated B~V colors and Ha equivalent widths 
(EW (Ha)) using the model's known SFRs and mean oxy- 
gen abundances. The interested reader can find details of 
the population synthesis model and essential tests in the 
Appendix. There is observational evidence that LSB galax- 
ies are characterized by B — V colors and EW(Ha) that lie 
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in certain limits. For instance, the observed B — V colors 
of blue LSB galaxies (which we attempt to model in this 
pape r) appear to lie in the 0.3 — 0.7 range (de Blok et al.l . 
119951 ). According to IZackrisson et all ((2005j), EW(Ra) of 
a sample of blue LSB galaxies seem to be confined within 
rather narrow bounds, EW(Ra) = 10-60 A. Hence, popu- 
lation synthesis modelling may be regarded as a consistency 
test for our numerical models, which should yield EW(iia) 
and B — V colors in acceptable agreement with observa- 
tions. 

We present the results of population synthesis mod- 
elling only for models with different IMFs; models with 
different SFRs show a similar behaviour. Figure [TU] shows 
model EW(Ha) (left column) and B — V colors (right col- 
umn) as a function of galactic age for model 1 (top) and 
model 2 (bottom). The equivalent widths are derived us- 
ing SFRs averaged over the past 20 Myr. For convenience, 
the thick dashed lines show EW(Ha) averaged over 1 Gyr. 
The horizontal dotted lines in the left pa nels show a typical 
range for EW(Ra) in blue LSB galaxies (jZackrisson et all 
120051 ). while the horizontal dotted lines in the right pan- 
els outline typ ical B — V colors of blue LSB galaxies 
(|de Blok et all 11995ft . The arrow shows a (possible) red- 
dening due to dust extinctio n for a fiducial value of the 
color excess E(B -V) = 0.1 (|McGaughl . [l994h . 

It is evident that model EW(Ha) show considerable 
fluctuations of order 15-20 A around mean values, which 
reflect the corresponding fluctuations in the SFR. Nether 
model 1 nor model 2 can account for the observed equiv- 
alent widths in the early evolution, yielding too large 
EW(Ha) during the first 5-6 Gyr. Both models begin to 
match the observed EW (Ha) only in the late evolution af- 
ter 6 Gyr, with model 2 (Salpeter IMF) yielding a somewhat 
better fit. We note that all models feature short-term varia- 
tions in EW(Ha) of such an amplitude that the correspond- 
ing equivalent widths may accidently go beyond either the 
upper or lower observed limits in the late evolution phase. 
This is however not alarming, since these fluctuations are 
short-lived and the observed range is inferred from just a 
few galaxies and may enlarge as more observational data 
become available. 

Unlike model EW(Ha), there are no large short-term 
fluctuations in model B — V colors shown in the right col- 
umn of Fig. 1101 Younger galaxies have on average bluer 
colors. The lowest observed color B — V = 0.3 is reached 
after approximately 2.5-3.0 Gyr of evolution. Taking into 
account a possible reddening, the latter value may decrease 
to approximately 1.5-2.0 Gyr. A similar lower limit on the 
galactic age was inferred from numerical modelling of oxy- 
gen abundances in Section UTTl However, model equivalent 
widths suggest a much larger value, ^5 — 6 Gyr. This 
value may decrease slightly, if LSB galaxies feature trun- 
cated IMFs with a smaller upper mass limit, say 40 M 
instead of 100 M Q . The oxygen abundance fluctuation spec- 
trum T{A) can be used to better constrain the minimum 
age of blue LSB galaxies. 

It is well known that blue LSB galaxies have B — V col- 
ors that are systematica lly bluer than those of HSB galax- 
ies (|de Blok et al.l . 119951) . It is tempting to assume that this 
fact may be indicative of blue LSB galaxies being systemat- 
ically younger than their HSB counterparts. However, our 
evolved galaxies seem to have colors that fall within the 
observed limits. There is a mild possibility for the model 
B — V color to exceed the upper observed limit if dust ex- 
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Fig. 11. Mean oxygen abundance versus Ha equivalent 
width diagram. Each symbol presents a pair of data points 
corresponding to a model galaxy of specific age, as indi- 
cated in Gyr. Each symbol type corresponds to a partic- 
ular numerical model specified in the upper-left corner. 
Bidirectional bars assigned to every symbol characterize a 
typical range of values found in our numerical modelling. A 
gray rectangular area delineates a range of observed values 
typical for blue LSB galaxies. 



tinction is taken into account. But considering the uncer- 
tainty in dust content in blue LSB galaxies, this possibility 
is not conclusive. In any way, even if blue LSB galaxies are 
younger than their HSB counterparts, the difference in age 
is not dramatic, of the order of a few Gyr at most. It is likely 
that blue LSB galaxies owe their bluer than usu al colors to 
the sporadic nature of star formation (see e.g. iMcGaughl . 
11994 Ide Blok et all 119951 : iGerretsen fc de Blokl . ll999l) . 



6.3. Mean oxygen abundances versus Ha equivalent widths 

As the final stroke of the brush, we present the ([O/H]) — 
£W(Ha) diagram in Fig.QT] Each symbol in this diagram 
represents a pair of data points [([O/H]), £W(Ha)] corre- 
sponding to a specific galactic age in Gyr as indicated in the 
figure. In particular, filled circles and open squares repre- 
sent models 1 and 2, respectively. Each symbol is assigned 
bi-directional bars, which characterize a typical range of 
values found in our numerical simulations. More specifi- 
cally, the vertical bars illustrate typical short-term fluctu- 
ations of EW(Ha) within each model, while the horizontal 
bars represent a typical range of ([O/H]) found among mod- 
els with identical IMFs but different SFRs. A gray rectan- 
gular area shows a typical range of observed oxygen abun- 
dances and Ha equivalent widths in blue LSB galaxies. 

There are several important features in Fig. QT] that we 
want to emphasize. First, only those model galaxies that 
arc older than 5-7 Gyr appear to fit into the observed 
range of both EW(Ha) and B — V. Younger galaxies fail 
to reproduce the observed Ha equivalent widths. Second, 
numerical models with standard mass limits (mi ow ,m up = 
0.1,100 Mq) populate the upper-right portion of the ob- 
served EW(Ha) vs. B — V phase space. This suggest that 
at least some LSB galaxies may have IMFs that are different 
from those considered in the present work, i.e., they may 
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have a truncated upper mass limit. Third, our modelling 
suggests an upper limit on the age of blue LSB galaxies 
that is comparable to or slightly less than the Hubble time. 
Forth, an accurate age prediction is possible if some inde- 
pendent knowledge of the IMF is available. 



7. Conclusions 

We have used numerical hydrodynamic modelling to study 
the long-term (~ 13 Gyr) dynamical, chemical, and pho- 
tometric evolution of blue LSB galaxies. Motivated by ob- 
servational evidence, we have adopted a sporadic model for 
star formation, which utilizes highly localized, randomly 
distributed bursts of star formation that inject metals into 
the interstellar medium. We have considered several model 
galaxies that have masses and rotation curves typical for 
blue LSB galaxies and are characterized by distinct SFRs 
and IMFs. We seek to find chemical and photometric signa- 
tures (such as radial oxygen abundance fluctuations, time 
evolution of mean oxygen abundances, integrated B—V col- 
ors and Ha equivalent widths) that can help to constrain 
the age of blue LSB galaxies. We find the following. 

— Local bursts of star formation leave signatures in the 
disk in the form of oxygen enriched patches, which are 
subject to a continuing destructive influence of differ- 
ential rotation, tidal forcing from spiral stellar density 
waves, and supernova shock waves from neighboring star 
formation sites. As a result, our model galaxy reveals 
short scale (~ 1-2 kpc) radial fluctuations in the oxygen 
abundance [O/H]. Typical amplitudes of these fluctua- 
tions are 0.5-1.0 dex in (l-2)-Gyr-old galaxies but they 
decrease quickly to just 0.1-0.2 dex in (5-13)-Gyr-old 
ones. 

— The oxygen abundance fluctuation spectrum (i.e., the 
normalized number of fluctuations versus the fluctua- 
tion amplitude) has a characteristic shallow slope in 
(l-2)-Gyr-old galaxies. The fluctuation spectrum is age- 
dependent and steepens with time, reflecting the ongo- 
ing mixing of heavy elements in the disk. The depen- 
dence of the fluctuation spectrum on the model param- 
eters and spatial resolution makes it possible to use the 
spectrum only for order-of-magnitude age estimates. 

— The mean oxygen abundance ([O/H]) versus EW(Ha) 
diagram can be used to constrain the age of blue LSB 
galaxies, if some independent knowledge of the IMF is 
available. 

— Numerical models with Salpeter and Kroupa IMFs and 
standard mass limits (0.1-100 Mq) populate the upper- 
right portion of the observed ([O/H]) vs. EW(Ha) 
phase space. This implies that blue LSB galaxies are 
likely to have IMFs with a truncated upper mass limit. 

— Our model data strongly suggest the existence of a min- 
imum age for blue LSB galaxies. From model B — V 
colors and mean oxygen abundances we infer a mini- 
mum age of 1.5-3.0 Gyr. However, model Ha equivalent 
widths imply a larger value, of order 5-6 Gyr. The lat- 
ter value may decrease somewhat if LSB galaxies host 
IMFs with a truncated upper mass limit. A failure to 
observationally detect large oxygen abundance fluctua- 
tions, which, according to our modelling, arc typical for 
(l-2)-Gyr-old galaxies, will argue in favour of the more 
evolved nature of blue LSB galaxies. 



— The oldest galaxies (13 Gyr) in all numerical models 
have mean oxygen abundances, colors and equivalent 
widths that lie within the observed range, suggesting 
that LSB galaxies are not considerably younger than 
their high surface brightness counterparts. 

— Low metallicities and moderate cooling render the for- 
mation of cold gas with temperature Sa 50 K inefficient. 
This fact implies that most of the cold phase may be 
in the form of atomic rather than molecular hydro- 
gen, though further numerical simulations are needed 
to make firm conclusions. 

— Sporadic star formation yields no radial abundance gra- 
dients in the disk of our model galaxy, in agreement 
with the lack of oxygen abundance gradients in blue 
LSB galaxies. 
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Appendix A: Disk scale height 

We derive the disk vertical scale height Z at each computational cell 
via the equation of local vertical pressure balance 



Pg c s 



dz, 



(A.l) 



where p g is the gas volume density, gf a " and gj* are the vertical grav- 
itational accelerations due to self-gravity of gas and stellar layers, 
respectively, and g!? al0 is the gravitational pull of a dark matter halo. 
Assuming that p g is a slowly varying function of vertical distance z 
between z = (midplane) and z = Z (i.e. E g = 2 Z p g ) and using the 
Gauss theorem, one can show that 
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where r is the radial distance and M^{r) is the mass of the dark matter 
halo located inside radius r. In equation l lA.31 we have assumed that 
the volume density of stars can be written as S s t/Z. Substituting 
equations flA.21 -l lA. 11 back into equation JA. Il l we obtain 
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This can be solved for p g given the model's known c^, S g , S s t, and 
Afj 1 (r). The vertical scale height is finally derived as Z = £ g /(2p g ). 



Appendix B: Population synthesis modelling 



We make use of the single-burst models by Bruzual & Chariot (2003) 
to synthesize broad-band colors and Ha equivalent widths. The con- 
tinuum level F c (t) near the Ha line at time t can be found as 



Fc(t) 



f 



f c (t-T,Z m )H>(T)dT 



(B.l) 
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Fig. B.l. Ha equivalent width versus galactic age for a 
test model as described in the text. The solid line presents 
the data for metallicity Z m = 0.001. The horizontal dotted 
lines delineate the observed range of values for blue LSB 
galaxies. 



where ^(t) is the SFR at time t, and f c (t — t, Z m ) is the continuum 
luminosity per wavelength for a single-burst stellar population with 
age t — t and metallicity Z In . 

The Ha luminosit y Ln a (t) at time t is estimated as Ljj a [erg/s] = 
1.36 X 10~ 12 JV Lyc [l/s] (IKennicuttl . 119831) . where iV Lyc is the rate of ion- 
izing photon production derived directly from Bruzual and Chariot's 
population synthesis modelling. The resulting equivalent width is 
found as EW(Ha) = Ln a /F c . The broad-band B — V color at time 
t is found via B and V luminosities (Lb(*) and Ly(t), respectively) 
calculated using equation J B .11 ). in which f c is replaced with lumi- 
nosities in the B and V bands, respectively. The resulting color is 
(B-V) = -2.51og 10 [L B (t)/L v W]. 

We perform a consistency test on our populatio n synthesis model 
by com paring our predicted EW(Ha) with those of FZackrisson et al.l 
(2005) for a model galaxy with a constant (over 15 Gyr) SFR, total 
mass of formed stars 10 10 Mq, metallicity Z m = 0.001, and Salpeter 
IMF with mass limits (rrii ow , m up ) = (0.1, 120) M®. The resulting Ha 
equivalent width as a function of time is shown in Fig.I B.l l bv the solid 
line. O ur model Ha equivalent widths match those of Zackrisso n et al.l 
(120051 ) quit e well— at t = 15 Gyr w e find EW(Ha) = 88 A, whereas 
figure 2 in IZackrisson et al.l (120050 shows EW(Ha) « 86 A. This 
example demonstrates that a constant SFR produces too large Ha 
equivalent widths and hence is unlikely. 

Finally, we would like to note that Z m = 0.001 (or [O/H] = 
— 1.3) is not typical for p r esent- day b lue LSB galaxies. Acco r ding t o 
Ide Blok fc van der HuTstl ( 119981) and IKuzio de Narrav et all (120041) , 
the mean oxygen abundance of blue LSB galaxies is ([0/H])lsb = 
—0.85 dex or one seventh of the solar metallicity. Higher metallicity 
is expected to yield lower Ha equivalent widths. 
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